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Abstract 

We propose a new law for the deceleration parameter that varies linearly with time 
and covers Berman's law where it is constant. Our law not only allows one to generalize 
many exact solutions that were obtained assuming constant deceleration parameter, 
but also gives a better fit with data (from SNIa, BAO and CMB), particularly con- 
cerning the late time behavior of the universe. According to our law only the spatially 
closed and flat universes are allowed; in both cases the cosmological fluid we obtain 
exhibits quintom like behavior and the universe ends with a big-rip. This is a result 
consistent with recent cosmological observations. 
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1 Introduction 

It has been announced in 1998 that type la Supernova observations indicate that the current 
universe is not only expanding but also accelerating. This was perhaps the most striking 
discovery of the modern cosmology. Today, this behavior of the universe is confirmed by 
various independent observational data, including the Type la Supernova (SNIa), the large 
scale structure, the cosmic microwave background (CMB) radiation, and so on. There is a 
consensus on the conclusion that the universe has entered a state of accelerating expansion 
at redshift z ~ 0.5. In two recent studies by Cunha and Lima [1] and Cunha [2j, the 
transition redshift to the accelerating expansion of the current universe has been examined 
using the kinematic approach to cosmological data analysis that provides a direct evidence 
to the present accelerating stage of the universe which does not depend neither on the 
validity of general relativity, nor on the matter-energy content of the universe. Today 
we can examine not only when the cosmic acceleration began and the current value of 
the deceleration parameter, but also how the acceleration (the deceleration parameter) 
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varies with time. In a recent paper, Li et al. ^ has examined the current acceleration of 
the universe by the largest and latest SNIa sample (Union2), baryonic acoustic oscillation 
(BAO) and cosmic microwave background (CMB) radiation together. They conclude that 
once the systematic error is taken into account, two different subsamples of Union2 along 
with BAO and CMB all favor an increase of the present cosmic acceleration. The authors 
also give plots they obtained for the change of the deceleration parameter q with redshift 



Motivated from the studies outlined above, we propose here a linearly varying decelera- 
tion parameter (LVDP), which can be used in obtaining accelerating cosmological solutions. 
As a special case, LVDP also covers the special law of variation for Hubble parameter, which 
yields constant deceleration parameter (CDP) models of the universe, presented by Berman 
[H [5] . LVDP is consistent with the results of Cunha, and exhibits a similar behavior to 
the one obtained by Li et al. [3] when we plot the deceleration parameter as a function of 
redshift. After the discovery of the late time acceleration of the universe, many authors 
have used CDP to obtain cosmological models in the context of dark energy (DE) in general 
relativity and some other modified theories of gravitation such as f{R) theory within the 
framework of spatially isotropic and anisotropic spacetimes. However, generalizing CDP 
assumption would allow us to construct more precise cosmological models. Indeed, as we 
show in this paper, a simple generalization of CDP ansatz to LVDP ansatz fits the recent 
observational results pretty well. In this paper, for convenience, we consider only the spa- 
tially homogenous and isotropic Robertson- Walker (RW) spacetimes and hence the perfect 
fluid representation for the energy-momentum tensor in the context of general relativity. 
The LVDP ansatz can also be used in spatially homogenous but anisotropic spacetimes 
and hence the imperfect fluids might also be considered in the context of general relativity, 
as well as generalized theories of gravity such as f{R) gravity. 

2 The field equations, law for the deceleration parameter 
and solution 

The Einstein field equations can be written as follows. 



where G^u is the Einstein tensor and T^^ is the energy-momentum tensor. To solve the 
Einstein field equations it is usually necessary to make some simplifying assumptions such 
as choosing a metric with a significant degree of symmetry. In this paper, we consider the 
Robertson- Walker metric with a maximally symmetric spatial section. 
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where a{t) is the cosmic scale factor and the spatial curvature index k = —1, 0, 1 corre- 
sponds to spatially open, flat and closed universes, respectively. Considering a co-moving 
fluid, the RW metric admits only the perfect fluid representation for the energy-momentum 
tensor, which can be written as follows, 

Ti: = [p,p,p,p] (3) 

where p = p{t) is the energy density and p = p{t) is the pressure of the fluid. In a comoving 
coordinate system, the Einstein field equations (1) with (2) and (3), read as 

3- + 3- = p, 4 

a} a a 

— + 2- + — = -p. 5 
a a'^ 

This system consists of two equations but three unknowns (a, p and p), hence is not fully 
determined. One may determine the system fully by specifying a theory that determines 
a relation between the energy density and pressure of the fluid. Most of the perfect fluids 
relevant to cosmology obey an equation of state of the form 

p = wp, (6) 

where w is the equation of state (EoS) parameter, not necessarily constant. The three most 
common examples of cosmological fluids with constant w are the dust {w = 0), radiation 
{w = ^) and vacuum energy {w = —1) that is mathematically equivalent to cosmological 
constant A. It is well known that fluids with w < are usually considered in the context 
of DE, since they give rise to accelerating expansion. Besides the fluids with constant 
EoS parameter with w < particularly in the context of DE, there have been proposed 
various scalar fleld models that can be described by a time dependent w that can evolve 
below the — |, e.g., quintessence (— 1 < w < 1), phantom {w < —1), quintom that can 
evolve across the cosmological constant boundary w = —1. However, we still do not have 
a perfect understanding of the EoS of the DE and sometimes it is useful to think about 
the Einstein fleld equations without specifying the theory from which EoS parameter w is 
derived [6]. Hence, as an alternative, we may seek for the solutions of the Einstein fleld 
equations under an ansatz or some ansatze concerning the kinematics consistently with 
the observed kinematics of the universe. Then one can derive an EoS for the possible 
fluid that might be driving the expansion of the universe and particularly the current 
acceleration of the universe in the context of DE. We may, then, check whether such a 
fluid is a realistic energy-momentum source or not and try to learn about the nature of 
the DE. Since the twice contracted Bianchi identity {G^'^-u = 0) assures the conservation 
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of the energy-momentum tensor (T^'^-i, = 0), we would be sure that the fluid would satisfy 
the energy-momentum conservation law. Furthermore, we should require positivity of the 
energy density of the co-moving fluid p > to be satisfied as a condition for our fluid to 
be a realistic energy-momentum source. 

Following a similar line of reasoning, Berman [1] and Berman and Gomide |5j proposed 
in the 1980's a law of variation for Hubble parameter within the context of RW spacetimes 
in general relativity that yields constant deceleration parameter {q = = m — 1, where a 
is the scale factor and m > is a constant). In Berman's law the deceleration parameter can 
get values q > —I, and since —l<q<0 corresponds to the accelerating expansion, many 
authors have studied cosmological models using this law in the context of DE following the 
discovery of current acceleration of the universe. In this paper, we propose a generalized, 
linearly varying deceleration parameter 

aa 

q = -— = -kt + m - 1, (7) 

where A: > and m > are constants, k = reduces ([7]) to the law of Berman [H [5]. 
Using this law one can generalize the cosmological solutions that are obtained via CDP. 
Furthermore, one can obtain models that fit better the cosmological data. The universe 
would exhibit decelerating expansion if g > 0, an expansion with constant rate if g = 0, 
accelerating power-law expansion if — 1 < q < 0, exponential expansion (also known as de 
Sitter expansion) if g = —1 and super-exponential expansion if g < — 1 (see fSjlUj for super- 
exponential expansion). Hence, the fastest rate of expansion under the CDP ansatz is the 
exponential expansion. However, with our LVDP ansatz the universe inevitably evolves 
into the super-exponential expansion phase unless k = 0, which is one of the possible fates 
of the universe according to the cosmological observations (see |6| UH [12]). 
Solving ([7]) one obtains three different form of solutions for the scale factor: 



for /c> and m > 0, (8) 

a = 02 (mt + 02)™ for A; = and m > 0, (9) 

a = ase"^^* for A; = and m = 0, (10) 

where oi, 02, 03, ci, C2 and C3 are constants of integration. The last two of these solutions 
are for constant q, and hence correspond to the solutions under CDP ansatz. We will not 
dwell on these but only on the first one, which is new. For convenience, in the following we 
consider the solution for A; > and m > and omit the integration constant ci by setting 
ci = 0. By doing this, we also set the initial time of the universe to ti = 0. The reason 
for considering the solution only for A; > and m > is not only for simplicity but also 
for compatibility with the observed universe. A; > means we are dealing with increasing 
acceleration {q = —k < 0). Because tj = and A; > 0, the only way to shift the deceleration 
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parameter to values higher than —1 is to set m > 0. We would like to note that here, before 
proceeding, if one would like the universe to commence with a decelerating expansion, one 
may tighten our ansatz by choosing m > 1. Under the above considerations, ^ is further 
reduced to 

„ — arctanhf — t-l) /-i 1 ^ 

Now, using this, the Hubble parameter of the universe is obtained as follows: 

H = ^ = - „ ' (12) 

a t{kt - 2m) ^ ' 

Substituting the scale factor pT]) in dH and ^ we get the energy density 
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and the pressure 

^ t^{kt-2my ai^ ^ ' 

Finally, using (fT3l) and ([H]) in ([6]) the EoS parameter of the fluid is obtained as follows: 



-8ai2(A;t-m + |)e^'^=-^*'^°K^*-i) -Kt2(A:t-2m)2 ^ ^ 

W = T TT T . (15) 

12ai2 e™'^''"*^"H^*-i) + SKf^ikt - 2m)2 



We also solve for the deceleration parameter q asa function of the redshift z = — 1+^^^, 
where 02=0 is the present value of the scale factor: 



j{z) = 2m — 1 — mtanh 



— ln(z + 1) — arctanh ( — ~^ ^ _ 2 
2 ^ ^ \ m 



(16) 



We know that qz = qz=o = — 1 + m for CDP (i.e., the A: = case in our model). Indeed, 
one may check that the choice k = reduces (|16p to q{z) = —1 + m. Note that the above 
relation is independent of the spatial curvature. 



3 Physical behavior of the model and observational valida- 
tion 

In our model the universe has finite lifetime. It starts with a big bang at ti = and ends 
at iend = Both of the energy density of the fluid and the scale factor diverge in finite 
time as t — )• tend- This is the big rip behavior first suggested by Caldwell et al. [11]. The 
universe begins with = m — 1, enters into the accelerating phase {q < 0) at tt = ^^^^^ 
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(assuming m > 1), enters into super-exponential expansion phase {q < —1) at tge = ^ and 
ends with ggnd = — 1- 

One may observe that positivity condition on the energy density of the fluid is valid 
for the spatially closed and flat universes but violated for the spatially open universe (the 
reason being that k = — 1 < 0). Hence, under the ansatz of LVDP ([7]), spatially closed and 
flat universes are possible while spatially open universe is impossible. It may be worthwhile 
to note that Berman [1] also reached the same conclusion. 

The EoS parameter of the fluid exhibits different behaviors for spatially closed, flat and 
open models. However, in all of the three models w = at t = tt {as would be expected) 
and w = —\rn — 1 at t = tend- 

Within the context of general relativity, the current cosmological data from SNIa (Su- 
pernova Legacy Survey [13], Gold sample of Hubble Space Telescope [H]), CMB (WMAP, 
BOOMERANG) [K] and large scale structure (SDSS) [IB] do not rule out DE with a dy- 
namical equation of state parameter that can evolve into the phantom region [w < —1). El 
This is consistent with our model. Additionally, because the observational data allows DE 
to pass into the phantom region, it means that the data also do not rule out the possibility 
of a big rip in future universe, which is also consistent with our model. 

From model dependent or independent analyses of the cosmological observations in the 
literature [H [21 [Sj [U [171 (iBl IHl ISO] , the transition redshift of the accelerating expansion 
is given by 0.3 < Zt < 0.8. In particular, the kinematic approach to cosmological data 
analysis provides a direct evidence to the present accelerating stage of the universe, which 
does not depend on the validity of general relativity, as well as on the matter-energy content 
of the universe [1]. By assuming a spatially flat universe Cunha [I] finds that qz=o = —0.73 
and Zt = 0.49. Now, to demonstrate how LVDP matches the observed kinematics of the 
universe and makes additional predictions, we first plot the cosmological parameters by 
choosing m = 1.6 and k = 0.097. 

In Fig. [Tjwe plot the scale factor a versus cosmic time t. The universe starts with big 
bang at t = and ends at tend — 33, which is close to the given lifetime for the universe 
by Caldwell |llj trip = 35. In Fig. [2] we plot the Hubble parameter H versus cosmic time 
t. The Hubble parameter diverges at the beginning and end of the universe. In Fig. [3] 
we plot the deceleration parameter q versus cosmic time t. The deceleration parameter is 
initially qi = 0.6 and reaches gend = —2.6 at the end of the universe. The universe enters 
into the accelerating expansion phase at t = 6.2 and the present (i.e. at t = 13.7) value of 
the deceleration parameter is q = —0.73. Both values are consistent with the observational 
results. 

One may observe that the spatial curvature parameter k is not included in the kine- 
matics but in the dynamics of the fluid. This is what we expect, because we constrain the 
kinematics from the beginning but not the dynamics of the fluid. Hence, we present plots 
of the energy density, pressure and EoS parameter of the fluid with different colors for the 

^See [6l [T] m [12] for theoretical and observational status of crossing the phantom divide line. 
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spatially closed (green), flat (red) and open (blue) universes. In Fig. [3]we plot the energy 
density of the fluid p versus cosmic time t. One may observe that the spatially open model 
is not possible since the positivity condition of the energy density is violated but the spa- 
tially closed and flat cases are possible since the positivity of the energy density is satisfied 
in these models. Both for the spatially closed and flat models the energy density of the fluid 
diverges at the beginning and end of the universe. However, they behave differently and 
the energy density of the fluid in the spatially closed model is always higher than the one 
in the flat case. Because both of the scale factor and the energy density of the fluid diverge 
at tend) we say that the universe ends with a big rip. In Fig. [5] we plot the pressure of the 
fluid p versus cosmic time t. The pressure also diverges at the beginning and end of the 
universe but exhibits different behaviours for spatially closed, flat and open universes. In 
Fig. [6] we plot the EoS parameter of the fluid w versus cosmic time t. It exhibits different 
behaviors for spatially closed, flat and open universes, but converges to the same value at 
the very late times of the universe, strictly speaking, w — )-~ —2 as t — 7>~ 33. The vertical 
line in the graph for the spatially open model is there since the energy density in spatially 
open model changes sign at that point. Hence, as already mentioned above, this solution 
is not possible since it violates positivity of the energy density. 

We also plot the deceleration parameter q versus redshift z in Fig. [7]with solid curve by 
using the same values we used in the other graphs. One may observe that the accelerating 
expansion begins at 2; = 0.5 consistent with the observational data. We particularly would 
like to demonstrate that by choosing m = 1.45 and k = 0.08 we can obtain the plot of 
q{z) (the dashed curve in Fig. [7]) which fits quite well to the plots of q{z) for spatially flat 
universe given by Li et al. [3] in a recent paper, where they analysed the data from SNIa, 
BAO and CMB simultaneously. The bottom line argument of these plots is that different 
data sets and analysis methods may give different values for qz=Q and tt but the general 
behavior of q{z) follows a similar pattern. The pattern same is obtained via the LVDP 
ansatz that fits the observational curves quite closely. 

4 Discussion and conclusion 

We propose a special law for the deceleration parameter that is linear in time with a negative 
slope. The law we suggest also covers the law of Berman [Jl [5] where the deceleration 
parameter is constant. Berman's law has been used for obtaining exact cosmological models 
in the literature, particularly in the context of dark energy following the discovery of the 
current acceleration of the universe. The law we suggest gives the opportunity to generalize 
many of these models with better consistency with the cosmological observations. The big- 
rip predicted in our models also is not ruled out by observations. The models are possible 
only for spatially closed and flat universe and in both models the fluid we obtained exhibits 
quintom like behavior. 

Berman's law has also been applied to the Bianchi I spacetime assuming that the 
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energy-momentum tensor itself is isotropic. Later on, it has been applied to various other 
spacetimes assuming isotropic energy-momentum tensor by many authors. Considering 
anisotropic spacetimes, one may also generalize the energy- momentum tensor in a way so 
as to yield anisotropic pressure. For instance, Akarsu and Kilinc [21j and Kumar and Singh 
|22] applied the law presented by Berman to Bianchi I spacetime by allowing anisotropy in 
the pressure of the fluid. Similarly, the law we present can also be used within the frame- 
work of spatially homogeneous but anisotropic Bianchi type spacetimes and Kantowski- 
Sachs spacetime in the presence of isotropic and/or anisotropic fluid. To do that, one 
may define the mean scale factor, most generally, as S = {abc)^ , where a, b and c are the 
directional scale factors, and then generalize ([7]) as follows, 

SS 

q = —^ = -kt + m-l, (17) 

where A; > and m > 0. 

The law we have proposed can lead to many applications in cosmology; not only in 
the context of general relativity but also in generalized gravity theories such as the f{R) 
theory of gravity. 
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Figure 1: Scale factor a versus cosmic time t. The curve is plotted by choosing m = 1.6 
and k = 0.097. a diverges as t — >-~ 33. 
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Figure 2: Hubble parameter H versus cosmic time t. The curve is plotted by choosing 
m = 1.6 and k = 0.097. 
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Figure 3: Deceleration parameter q versus cosmic time t by choosing m = 1.6 and k = 0.097. 
The values of m and k are choosen by considering the results of Cunha [2] that obtained 
by kinematic analysis. 
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Figure 4: The energy density of the fluid p versus cosmic time t for the spatially closed 
(green), flat (red) and open (blue) universes by choosing m = 1.6 and k = 0.097. 
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Figure 5: The pressure of the fluid p versus cosmic time t for the spatially closed (green), 
flat (red) and open (blue) universes by choosing m = 1.6 and k = 0.097. 
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Figure 6: The equation of state parameter w versus cosmic time t for the spatially closed 
(green), flat (red) and open (blue) universes by choosing m = 1.6 and k = 0.097. Dashed 
line is w = — |. The vertical line in the pl(2t for spatially open universe is because the 
energy density of the fluid changes sign at that time. 






Figure 7: Deceleration parameter q versus redshift z. The solid curve is plotted by choosing 
qz=o = —0.73 and m = 1.6 considering the kinematic analysis of Cunha [2]. The dashed 
curve is plotted by choosing gz=o = —0.65 and m = 1.45 considering the analysis of Li et 
al [3]. One may check that the dashed plot fit to the graphs given by Li et al. for spatially 
flat universe quite well. 
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